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In this review, we briefly overview recent works on hybrid (nano) opto-mechanical systems that contain both mechanical oscillators 
and diamond nitrogen-vacancy (NV) centers. We review two different types of mechanical oscillators. The first one is a clamped 
mechanical oscillator, such as a cantilever, with a fixed frequency. The second one is an optically trapped nano-diamond with 
a build-in nitrogen-vacancy center. By coupling mechanical resonators with electron spins, we can use the spins to control the 
motion of mechanical oscillators. For the first setup, we discuss two different coupling mechanisms which are magnetic coupling 
and strain induced coupling. We summarize their applications such as cooling the mechanical oscillator, generating entanglements 
between NV centers, squeezing spin ensembles and et al. For the second setup, we discuss how to generate quantum superposition 
states with magnetic coupling, and realize matter wave interferometer. We will also review its applications as ultra-sensitive mass 
spectrometer. Finally, we discuss new coupling mechanisms and applications of the held. 
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1 Introduction 

In 2012, David Wineland, together with Serge Haroche, won 
Nobel prize ’’for ground-breaking experimental methods that 
enable measuring and manipulation of individual quantum 
systems” m. In his lab, the motion and the electron spin of 
trapped ions couple with each other. Therefore, the motional 
degree of freedom of the trapped ions can be coherently ma¬ 
nipulated, such as the ground state cooling, Fock states gener¬ 
ating and detecting with unit fidelity and efficiency O. Using 
these technologies, David Wineland generated Shrodinger’s 
cat states with trapped ions ||3l|4l, which for many years only 
existed in through experiment. 
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The ’’true” Shrodinger’s cat state is defined by a micro¬ 
scopic two-level system entangling with a macroscopic sys¬ 
tem. New technology, such as opto-mechanics, is needed 
to generate quantum superpositions in macroscopic systems. 
Opto-mechanics studies the radiation pressure on the me¬ 
chanical oscillator, and how to use this pressure to manip¬ 
ulate the motion of mechanical systems. In 17th century Ke¬ 
pler found the radiation pressure of light when he studied the 
comet tails. In 1909, Einstein studied the statistics of radi¬ 
ation pressure force fluctuations on movable mirror O. In 
1960s, Braginsky studied the role of radiation pressure in the 
context of interferometers, and proposed to cool the motion- 
nal temperature of the mirror by radiation pressure lO. In 
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1969, A. Ashkin observed optical trapping of micron-sized 
particles in liquid Q. Two years later, optical levitation of 
glass spheres by an upward-propagating laser beam in both 
air and vacuum was demonstrated O . The optical tweezers 
now is widely used in atomic physics, chemistry and biology. 

Recent years, we witness great developments of opto- 
mechanics (Mil. We have achieved the quantum ground 
state of macroscopic mechanical oscillators by traditional 
cryogenic techniques in 2010 ifT^ . or by active cooling in 
2011 |[T3l[T4l. We have developed many technologies based 
on opto-mechanics. For example, we have transferred the sig¬ 
nals from light to mechanical oscillation, and vise versa ifTSl . 
We have realized optomechanically induced transparency 
m and generated squeezed lights in opto-mechanical sys¬ 
tems riTi . By using mechanical oscillator as an interface, 
quantum information has been able to exchange between mi¬ 
crowave superconducting circuits and optical lights dHii. 

In future we aim to achieve the similar controlling ability 
of macroscopic resonators as we did in the trapped ions. Then 
we can generate Schrodinger’s cat states with macroscopic 
objects ll2T]| . or even macroscopic living objects, such as virus 
Ea. As macroscopic quantum superposition states become 
bigger and bigger, we hope to identify the intrinsic decoher¬ 
ence mechanisms such as gravity induced decoherence ||23]| , 
Continuous Spontaneous Localization models II24II25L and et 
al.. For readers who are interested in the theory and develop¬ 
ments of optomechanics, please read these reviews ll9HTT]|. 

In order to generate the non-classical states in an opto¬ 
mechanical oscillator, strongly quadratic coupling between 
the oscillator and the optical mode is essential ll2T]| . How¬ 
ever, it is usually very challenging to get the large quadratic 
coupling in experiments. Inspired by the trapped ion experi¬ 
ments, the quadratic coupling can be replaced by coupling the 
mechanical oscillator to a two-level quantum system, which 
can be controlled well externally. In order to get the sta¬ 
ble macroscopic quantum states, the coherence time of both 
the mechanical oscillator and the spin 1/2 should be maxi¬ 
mized. There are several proposals that based on quantum 
dots ||26l[271, atoms ll28ll . superconducting qubits ifT^ and 
Nitrogen-vacancy centers (NV centers) Il29ti32l . Among all 
these proposals, the ones with NV centers are most attrac¬ 
tive. As NV centers have presented the coherence even at 
room temperature. 

NV centers in diamond are usually regarded as artificial 
atoms in solid systems. The diamond lattice consists of cova¬ 
lently bond carbon atoms, which makes diamond very stiff. 
The valence electrons in diamond have huge bandgap (5.48 
ev), which makes it transparent deep into the UV. In its lat¬ 
tice, nitrogen and vacancy are the most common defects. As 
shown in Fig. [T] an NV center consists of a nitrogen atom and 
a vacancy at the nearest neighboring site. The NV centers are 
usually negatively charged, possessing 6 electrons and spin 
5 = 1 in the ground state. 




Figure 1 Left figure shows lattice structure of diamond with nitrogen- 
vacancy center. Right one shows the negatively charged NV center electronic 
energy level structure. Electronic spin polarization and readout is performed 
by optical excitation at 532 nm and red fluorescence detection around 637 
nm. Ground-state spin manipulation is achieved by resonant microwave ex¬ 
citation. The ground-state triplet has a zero magnetic held splitting ^ 2.87 
GHz. Figure adopted from Ref. (21 • Copyright (2012) by Nature Publishing 
Group. 

NV centers are promising candidate system for quantum 
information processing. Single NV centers can be addressed 
using the confocal microscopy technique. The spin state of 
NV centers can be initialized and read-out at room temper¬ 
ature, respectively, by optical pumping and spin-dependent 
fluorescence. Furthermore, the spin state manipulation can 
be achieved by resonant microwave radiation. Because of 
the weak spin-orbital coupling of diamond material, and low 
concentration of the only spin-carrying isotope of carbon 
with natural abundance 1.1%, electron spin of NV centers in 
diamond has very long coherence time (~ ms) even at room 
temperature m. All these properties mentioned above make 
NV centers excellent candidate for quantum logic elements 
for quantum information processing. 

In addition to the applications in quantum information, NV 
centers are widely used as solid-state ultra sensitive magnetic 
field sensor. Again, because of the long coherence time, a tiny 
change of the magnetic field can be monitored by measuring 
its effect on the spin dynamics of NV centers. Recent exper¬ 
iments demonstrated high sensitivity ~ nT/ VHz at atomic 
scale resolution ||35l[36l. Furthermore, recent research shows 
that, not only magnetic fields, other types of signal (e.g., elec¬ 
tric fields, temperature, and strain etc.) which can be con¬ 
verted to magnetic signal can be also precisely measured by 
NV center spins. This makes NV centers very amazing multi¬ 
functional atomic-scale sensors. 

There are two methods to interface the NV centers and the 
mechanical oscillators. The first one is based on strong mag¬ 
netic field gradient to couple the mechanical oscillators with 
NV centers Ga. The second one requires the strain induced 
effective electric field to mix phonon mode with NV centers 
electron spins ll37]| . The strong coupling may reach in both 
of the methods. Once the strong coupling regime reaches, 
we may cooling the mechanical oscillator with coherent ex¬ 
citation exchange between spin and motion degrees. We can 
generate arbitrary superposition states of mechanical oscilla¬ 
tor. We may realize phonon laser and squeezing in the sys¬ 
tem |[371[38l. The system can also be used for quantum in- 
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formation processing (SOI, as it can interface many different 
degree of freedom. 

To increase the coherence time of mechanical oscilla¬ 
tor, we can optimize the design and choose the materials 
with high mechanical Q, such as Si 3 N 4 with Q around 10^ 
ll39l . However, the best way is optically trapping the nano¬ 
diamond in vacuum ll221l4Qfl48l . The mechanical Q factor in 
this system could be as high as 10^^, which is comparable to 
the trapped ions. It is possible to cool the trapped nanopar¬ 
ticle from room temperature B9l . Another advantage is that 
the strong coupling between the oscillation and the spin can 
be realized with a modest magnetic field gradient around 10^ 
T/m. It is also possible to couple the rotation degree of free¬ 
dom to the spins. Besides, the trapping frequency can be 
easily tuned, even completely turned off, which makes the 
time-of-flight measurement possible lEol. 

In an optical trap, the motion of nano-particle can be 
cooled to mK by feedback without cavity ||4T1[42. Ground 
state may reach by means of cavity sideband cooling II431I44L 
or by exchanging excitations between phonon and spin de¬ 
grees llisilollsll. Even at thermal states, matter wave inter¬ 
ference could be observed in this systems ||50l[52j[53l. By 
detecting the coherence of NV center spin, the ultra-sensitive 
mass spectrometer is realizable even at room temperature 

lEl. 

This review is organized as follows. In Sec. [2l we will 
review the methods and applications of strong coupling be¬ 
tween mechanical oscillator and NV centers. In Sec. [3l we 
will review how to observe macroscopic physics and achive 
ultra-sensitive detection in optically trapped nano-diamond 
that hosts NV centers. In Sec. IH we will summarize and 
discuss the potential developments in the field. 


2 Nanomechanical cantilevers coupled with 
NV centers 


To couple NV center electron spins with mechanical res¬ 
onator, we need to apply additional magnetic field gradient, 
or induce a strain in the lattice of diamond. Here we will 
briefly discuss these two different methods, and illustrate sev¬ 
eral applications for both of the methods. 


2.1 Magnetic field gradient induced coupling 

2.1.1 The setup 

zf \P . . 

' imagetictip 

^ i t 

NV center i 

I 

Figure 2 A magnetic tip, which is attached to the end of a cantilever, is 
placed at a distance h above a single NV center. It creats a strong magnetic 
gradient near the electronic spin of the NV center. As the spin energy is in 
proportional to the magnetic held, a strong coupling between the NV center 
electrons spin and the motion of cantilever is created. Laser held is used for 
initializing and measuring the spin states. Microwave is used for manipulat¬ 
ing the the spin states. 

The idea of strong magnetic coupling between a mechanical 
oscillator and a NV center is shown in Fig. I2lll29l. Here a 
single NV center electron spin is used for sensing the mo¬ 
tion of the mechanical oscillator. It should be noted that the 
NV center can also be attached to the end of resonator 1311 . 
other than the magnetic tip shown in Fig [2l The frequency 
of the mechanical oscillator is The oscillation induces a 
time-varying magnetic field which causes Zeeman shifts of 
the electron spin. The single phonon induces frequency shift 
X - gsPbGmCio, where ^ 2, yUg is Bohr magneton, Gm is the 
magnetic field gradient, and gq = AjhllmcOr is the zero field 
fluctuation for a resonator of mass m. We will show that in 
the current experiments, X can be as high as 100 kHz, which 
is larger than both the spin coherence time (1 ms), and the 
damping rate y = oJrIQ of high Q mechanical resonators. 

The Hamiltonian of the system can be written as 

Hs = H^v + hoJrO^a -f hX(a -f a^)Sz, (1) 

where H^v describes the dynamics of NV center elec¬ 
tron spin, a(a^) is the annihilation(creation) operator of the 
phonon in mechanical oscillator, and 5^ is the z component 
of the spin-1 operator of NV center. We assume that the NV 
sysmetry axis is also aligned along z direction. The NV center 
electron spin ground state is an 5 = 1 triplet. The three states 
are tabled as |0), | ± 1). Without external magnetic field, the 
levels I + 1) are degenerate. The zero field splitting between 
|0) and I + 1) is coqIIti = 2.88 GHz. In order to manipulate 
the levels | + 1) individually, we add a static magnetic field Bq 
along z axis to break the degenerate states | ± 1). The frequen¬ 
cies difference between |0) and | + 1) are 6 g)+i = (jjQ±pBBQlh. 

In order to coherent exchange the excitation between me¬ 
chanical resonator and the spin, we need to add microwaves 
to drive Rabi oscillations between |0) and | ± 1). We assume 
the magnetic field oscillation amplitudes is much less than 
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TlcoqIiib- In rotating wave frame we have 

Hnv = y + ^(|0>(/| + I0<0|), (2) 

where A+i and Q+i denote the detunings and the Rabi fre¬ 
quencies of the two microwave transitions. For simplicity, 
we assume that A+i = A_i = A and Q+i = = Q. 

The eigenstates of are |J) = (| - 1) - | -r 1))/V2, 
\g} = cos^lO) - sin^lZ?), and \e} = cos 6\b} sin^|0), where 

= (I + 1) + I - and tan(26i) = - V2Q/A. The 

eigenvalues corresponding to are = -A, and oje^g = 
(-A + VA^ -I- 2Q?)/2. In the case that A < 0, the lowest en¬ 
ergy state is |g). We adjust the frequency difference between 
\d) and \g) codg = oJd - ojg is equal to the mechanical fre¬ 
quency ojr, and the frequency difference between \e} and \d} 
<^ed = COe - Odd IS largely detuned with ojr- Under the condi¬ 
tion \(jJr - ojedl » T, we will have effective Jaynes-Cummings 
Hamiltonian 

Hjc = nAg\g}{e\aU h.c., ( 3 ) 

where Ag = -A sin 6. 

In this system, phonon-spin coupling strength A is usu¬ 
ally much less than phonon frequency cOy, which is around 
10 MHz. Therefore, it is impossible to directly generate 
Schrodinger’s cat states in this setup. We will discuss how 
to realize frequency tunable mechanical oscillator in the next 
section. By tuning the mechanical frequency less than A, 
the Schorodinger’s cat states in macroscopic mechanical res¬ 
onator can be prepared with spin-dependent force ll50l . 


2.7.2 Applications as resonator cooling 

With JC Hamiltonian Eq. (|3]), we can coherently manipulate 
the mechanical states with the help of spin states if strong 
coupling condition is fulfilled. The decay of mechanical res¬ 
onator at temperature T is defined as T^ = ksT I HQ. The 
spin dephasing rate is Scodg. The strong coupling condition 
is Ag > r^, Scodg. In experiments, we may use Si cantilever 
of dimensions (l,w,t) = (3,0.05,0.05) pm with frequency 
cOr = 1 MHz and = 5 x 10“^^ m. A magnetic tip can in¬ 
duce gradient field Gm = 10^ T/m at the distant h = 25 nm. 
The phonon-spin coupling strength T = 100 kHz is realizable. 
The mechanical Q can be as high as 10^ and the heating rate 
is r^/27r = 20 kHz. The spin decay can be around ScOgd = 1 
kHz if we use Carbon- 13 purified diamond. Therefore, the 
strong coupling condition can be fulfilled. 



Figure 3 Steady state phonon occupation number n/ as a function of en¬ 
vironment thermal excitation Nth- The value of denotes the crossover 
from LD regime to a regime where phonon cooling is strongly suppressed 
and finally rif = Nth- Firgure adopted from Ref. (21 • Copyright (2010) by 
The American Physical Society. 


To manipulate the phonon states of the resonator, the first 
step is to cool it to the ground state. If the mechanical os¬ 
cillator frequency is high enough, e.g. > 1 GHz, we can use 
the traditional cryogenic techniques to cool the environment 
to 10 mK, and the resonator is already in its quantum ground 
state. Usually the resonator frequency is around MHz, we 
have to use active method to cool it to ground state. In opto- 
mechanics, we need cavity mode as an vacuum bath to cool 
the mechanical mode. In this setup, we don’t need cavity, as 
NV center can be used as an effective vacuum bath. By using 
JC Hamiltonian (|3]) and optical pump technology of NV cen¬ 
ter spin states, we can cool the resonator to the ground state. 
The basic idea is as follows. We initialize the NV center spin 
to |0) by optical pump in the time scale less than 1 ps. Then 
we use microwave to prepare the spin state to |g). By turn¬ 
ing on the microwave driving we use JC Hamiltonian ^ to 
exchange transfer excitation from resonator to the NV center 
spin in the time scale l/Ag ~ 10 ps. As last we initialize the 
NV center again. By repeating the process, the excitations in 
resonator are quickly removed. Under the condition > T^, 
the quantum ground state cooling should be possible. If we 
carefully adjust the detuning of driving microwave and laser, 
the above procedure can automatically repeated to cool the 
mechanical resonator l|29l . 

The above cooling method is only valid for low tempera¬ 
ture r < 1 K and mechanical frequency ~ 1 MHz. For 
high temperature, the Lamb-Dicke(LD) limit is not fulfilled, 
the cooling with NV center spins is no longer possible Ea, 
as shown in Fig. [3l When thermal phonon phonon is larger 
than = 10^, the cooling effects quickly disappear. After 
quantum ground state being prepared, arbitrary superposition 
states of resonator should be possible. However, to generate 
the high fidelity phonon states, we need a high Q resonator, 
which is very difficult for this system. We will discuss this 
problem in Sec. [3] 
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2.1.3 Applications in quantum information processing 

Once mechanical ground state is prepared, the mechanical 
resonator can be used for quantum information processing 
ISbHSHl . By charging the mechanical oscillator, distant cou¬ 
pling between resonators is also possible 1301. 



Figure 4 A scheme to entangle multiple NV centers spins that couple with 
the same mechanical oscillator, (a) An array of evenly spaced magnetic tips 
are attached at the end of a nanomechanical cantilever, under which are one- 
to-one correspondent NV centers. Figure adapted from Ref. ED Copyright 
(2010) by The American Physical Society. 

In order to entangle distant NV centers, here we need to 
couple multiple NV centers to the same mechanical oscilla¬ 
tor, which is used as a quantum bus. As shown in Fig. [4^, Xu 
and et al. considered the scheme that an array of a NV centers 
locate equally distance, above which are one-to-one magnetic 
tips attached at the end a cantilever. The qubits are encoded 
in two spins levels 0) and -1). We applied a microwave to 
globally drive all the NV center spins. We assume that all 
NV center spins have the same transition frequency ojq and 
the frequency of microwave is on resonant with cjq. We de¬ 
fine the new basis for the spins as |+) = (| - 1) -f | -f 1))/ ^^2. 
Consider the case that the Rabi frequency of microwave Q is 
comparable to the mechanical frequency ojr, in rotating wave 
approximation, we have the effective Hamiltonian of the sys¬ 
tem as 

Hi = Y^ ^a\+)j{-\ exp[-/<5(fVf] + h.c., (4) 

with d{f) = ojr - n(0- The collective Hamiltonian Eq. dU) 
can be used for generating W states. As an example, initially 
we prepare the phonon state to the Fock state 11) with phonon 
number 1, and the NV centers to the state |-)i |-)2 • • • |-)n- 
We set the detuning S(t) = 0, the W state \W) = 

(l+)l|-)2-••!->. + |-)ll + )2---|->.---|-)l|-)2---|+).)/V^ 

can be generated at the time t = nl{ sfnA). 

If the detuning S(t) is much larger than coupling strength 
T, it will induce the effective spin-spin coupling between NV 
centers ||56l[59l. In Ref. (561, they only considered two NV 
centers coupling with the same cantilever. Assume the detun¬ 
ing Si = S 2 = S, to the second order we have the effective 
Hamiltonian 

Ml 1 

Heff= —[(a'a+-)(a-u + cr2z)-(a-i+a-2-+a-i-cr2+)]. (5) 

As we discussed in Sec. 12.1.21 the phonon-spin coupling HA 
could be around 100 kHz. In order to fulfill the condition 


d » Tg, we choose fiS in the order of 1 MHz. The spin-spin 
coupling fiA^I(4S) is in the order of 10 kHz, which is larger 
than both NV center spin decoherence rate and the effective 
phonon decay rate at 10 mK temperature. Therefore, it is pos¬ 
sible to generate the entanglement among NV centers spins in 
experiments with high fidelity. 



Figure 5 Schematic view of a scalable quantum information processor on 
an electromechanical quantum bus and NV centers spins as qubits. Figure 
adapted from Ref. Eq). Copyright (2010) by Nature Publishing Group. 

The idea can be extended to realize a scalable quantum 
information processor (SOj. As shown in Fig. [51 the proces¬ 
sor consist an array of N nanomechanical resonators, each of 
which is coupled magnetically to an electronic spin qubit as¬ 
sociated with an impurity located in the substrate below. As 
we discussed in Sec. 12.1.11 the resulting spin-phonon cou¬ 
pling could be as high as 50 kHz under the magnetic gradient 
10^ T/m. The long range interaction between distant sites 
can be established by charging the resonator and capacitively 
coupling with the nearby wire. The Hamiltonian of the res¬ 
onators is 

Her = X + I X = 2j 

i ij n 

( 6 ) 

After combining magnetic and electric coupling Hamiltoni¬ 
ans, we have the full Hamiltonian of the system. 

Her = Hs(t) + 7; X ^^nOlcin + ^ 2j (A 

n i,n 

where Hgt) = 2/ + 0^i(t)cr^f}l2, ojn and an denote fre¬ 

quencies and collective mode operators for phonon eigen- 
modes. This model is quite similar as the quantum com¬ 
putation proposal for trapped ions systems. However, here 
the decoherence mechanism and physical implementation are 
different. Anyway, the idea of using hybrid systems that con¬ 
tain mechanical resonators array and NV centers for quantum 
information processing needs to be further investigated in fu¬ 
ture. 
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2.1.4 Applications in ultrasensitive measurement 

As the resonator is coupled with NV center spin, we can 
use NV center as a detector to measure the state of res¬ 
onator EBEa. Therefore, the NV center spin can be used 
as an ultra-sensitive detector for mechanical resonators. Even 
strong coupling conditions is not fulfilled, we may still detect 
the Brownian motion of resonator in classical regime. 
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Figure 6 (a) A confocal microscope monitors the fluorescence of a single 

NV defect hosted in a diamond nanocrystal placed at the extremity of a SiC 
nanowire. A microwave antenna is used to manipulate the NV electronic 
spin, while a micro-fabricated magnetic structure approached in the vicin¬ 
ity of the suspended NV centre generates a strong magnetic held gradient. 

(b) , The electronic levels structure of the NV centres at zero magnetic held. 

(c) , Fluorescence map of the system recorded with the confocal microscope 
while scanning the objective position. The isolated bright spot circled in red 
corresponds to the fluorescence of a single NV centre. Inset: zoom on the 
nanowire extremity. Figure adapted from f3TI . Copyright (2011) by Nature 
Publishing Group. 


In Ref. ISH . it realized the hybrid system that a single 
nitrogen-vacancy center coupled with a nanomechanical os¬ 
cillator. A single NV center hosted in a nano-diamond is 
placed at the extremity of a Sic nanowire. By adding a strong 
magnetic field gradient to the system, a magnetic coupling 
between NV center and mechanical oscillator is induced. The 
motion of nanowire can be probed by reading out a single 
electron spin of the NV center. 

Reference El reported the coherent coupling of a me¬ 
chanical cantilever to the single spin of NV center. The au¬ 
thors used the electronic spin of an NV center to sense the 
motion of a magnetized cantilever. Although the magnetic 
field change due to the motion of the cantilever is very small, 
it can still be monitored by measuring the NV center spin co¬ 
herence. This was achieve with the help of the spin echo tech¬ 
nique, where the spin is Hipped by microwave n pulses during 
the evolution. With this technique, the unwanted background 
noise is filtered out, while the signal to be detected (i.e., the 
magnetic field signal from the cantilever motion) is magni¬ 
fied. In Ref. 1321 . the authors successfully demonstrated a 
motion sensing with precision down to ~ pm under ambient 
conditions. Later, we will show the spin-echo-based tech¬ 
nique is not limit in motion sensing. Indeed, it can be applied 
in sensing various types of signals lEl. 


2.2 Strain induced coupling 

2.2.1 Model 

In the previous subsection, we summarized the recent works 
on hybrid systems with a mechanical resonator and a NV cen¬ 
ter, where a strong magnetic gradient induces spin-phonon 
coupling. In this section, we will discuss new a mechanism 
that couples NV center electrons and the the phonon mode 
in diamond, strain induced spin-phonon coupling, and effec¬ 
tive spin-spin interactions between NV centers OTlEOl. The 
spin-spin interactions can induce spin squeezing in NV center 
ensembles (NVE), which may be used in NVE magnetome- 
try. 



Figure 7 (a) All-diamond doubly clamped mechanical resonator with an 

embedded NV centers ensemble (NVE). (b) Spin levels and transitions of NV 
centers electron spins. Local perpendicular strain induced by beam bending 
mixes the | + 1> states, (c) NV centers spins in the two-level subspace | + 1> 
off-resonantly couple with a common mechanical mode, (d) NVE collective 
spins squeezing. Figure adopted from Ref. 1^ . Copyright (2013) by The 
American Physical Society. 


As shown in Eig. [3a), we consider an NVE embedded 
in a single crystal diamond nanobeam. When beam oscil¬ 
lates, it strains the lattice of diamond, and induces direct cou¬ 
pling between NV center electrons spins and the mechanical 
mode. Here we write down the Hamiltonian for NV center, 
in the presence of external electric and magnetic fields E and 

i^NV =(^o + d\\Ez)Sl + pBgsS • B 

- dAEASxSy + SySd + E,(Sl - 5 ^)], 

where DqIIti = 2.88 GHz is the zero field splitting, gs - 2, 
yUfi is the Bohr magneton, and d\\ {dy_) is the ground state elec¬ 
tric dipole moment in the direction parallel (Perpendicular) 
to the NV axis (621. 

As shown in Eig. Ob), motion of nanoresonator changes 
the local strain of the NV center, and induces an effective 
electric field. We are interested to the near resonant coupling 
between a single motion mode and | + 1) transition of the NV 
center. The Zeeman splitting between | +1) is Ag = gsPsBzIli- 
The perpendicular component of strain E^ mixes | ± 1) states 
with interaction Hamiltonian E^ = E^^aA- a^), where a is the 
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destruction operator of the resonator mode with frequency 
(jjyn- Eo is the zero point motion of the beam induced per¬ 
pendicular strain. The parallel strain also induces the level 
shifts between | ± 1) and |0). However, in the subspace | ± 1) 
the parallel strain plays no role in the system, as state |0r) is 
not involved. In this subspace, the interaction Hamiltonian of 
each NV center is 


Hi = fig(a^cr- acr^), (9) 

where erf = | + l)(+l|is Pauli operator for ith NV center and g 
is the single phonon coupling strength. For NVE, we can de¬ 
fine the collective spin operators, ^ Z; |1)K1| “ I “ 
and J± = Jx± iJy = Z/ erf • The total Hamiltonian is 

Hp = CJmCL^CL AbJz + J- + ClJ+). (10) 

Here for simplicity, we assume that coupling g is uniform 
for all NV centers. The interaction between NV centers is 
also neglected as we assume that they are far apart from each 
other. 

The coupling strength g can be estimated with the follow 
equation El 


^ ^ 180( ^ _ ) 

2n 


1/2 


( 11 ) 


where p is the mass density and E is the Young’s modulus 
of diamond. Up to now, we only achieved the single phonon 
coupling strength g ~ 0.25 Hz in a diamond resonator m- 
Theoretically, we may get g ~ 1 kHz in future. For a beam 
of dimensions (L,w,h) = (1,0.1,0.1) pm we got 6L)^/27r ~ 1 
GHz and coupling g jin ~ 1 kHz. We define the single spin 
cooperativity rj = g^T 2 l(yhth), where y = is the me¬ 

chanical decay rate and 72 is the spin dephasing time, and fith 
is the thermal equilibrium phonon number for resonator. It is 
found that rf could be as high as 0.8 for 7 = 4 K, 72 = 10 ms, 
and Q = 10^. To get 7/ > 1, we need to further decrease the 
dimension of the diamond resonator, or decrease the temper¬ 
ature 7. Similarly as Sec., the Hamiltonian can also be used 
for cooling the mechanical oscillator if the strong coupling 
condition is fulfilled 13811 . 


2.2.2 Spin squeezing 

Similarly as we discussed in Sec. 12.1.31 |56l[59]l, in large de¬ 
tuning limit that g A = Ab - oJm, can adiabatically 
eliminate the mechanical mode a, and get effective spin-spin 
interactions El- The effective Hamiltonian is 

Hgff =■ (jjfYiCi^n (Ag -l- Aci^^Jz^ — 7+7_, (12) 

where A = 2g^/A. The Hamiltonian Eq. ([T^ can be used 
for generating spin squeezing state. To generate the spin 
squeezed state, we initialized the NVE to coherent spin state 
|j/^o) along X axis which satisfies 7jcl<Ao) = H has equal 

transverse variances (7^) = {J^} = Jjl. The interaction term 


J+J- = - 7^ -I- 7^, where total angular momentum 7 is 

conserved and 7^ induces spin variance squeezing in one di¬ 
rection, as shown in Fig. Hd). The squeezing can be qualified 
by the parameter 


.2 2/(A72.„) 

^ ’ 


(13) 


where (A7^-j^) = \(y+- + ^yz) is the minimum spin un¬ 

certainty with V± = {Jy ± Jj) and Vyz = {JyJz + JzJy)!^- If 
we prepare the spin squeezed state with < 1, it has a direct 
applications for magnetometry of NVE below the projection 
noise ilmit l6^ . 



Figure 8 (a) Spin squeezing parameter versus scaled precession time with 

= 100 spins. Thick blue (gray) lines show the calculated squeezing param¬ 
eter for T 2 = 10 ms and values of % as shown. For each curve, the squeezing 
is optimized. Dashed lines are calculated from the linearized equations for 
the spin operator averages, (b) Optimal squeezing versus number of spins. 
Lower (upper) line shows power law fit for iith = 1(10) and T 2 = 1(0.01) s 
. Other parameters in both plots are = 1 GHz, gl2n = 1 kHz, Q = 10^. 
Figure adopted from Ref. (^. Copyright (2013) by The American Physical 
Society. 


In Fig [It a), the squeezing parameter verse time is plot¬ 
ted for an ensemble of N = 100 spins and several fith- The 
decoherence of spin and mechanical decay is also consid¬ 
ered by solving master equation lITTl . We plot the scaling 
of the squeezing parameter with 7 for small decoherence in 
Fig. Ob). It was found the scalling is ~ 7“^^^. 


3 Optically trapped nano-diamond that hosts 
NV centers 


As we stated in the introduction, one of the key factors affect¬ 
ing macroscopic Schrodinger’s cat state is the mechanical Q 
factor. In this section we will focus on the optically trapped 
nanodiamond system, whose mechanical Q factor is not re¬ 
lated the the material, but the pressure of the vacuum. For the 
current technology, we can have Q has high as 10^^. There¬ 
fore, the life time of Schrodinger’s cat state can be as long as 
millisecond. 
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3.1 Scheme and Fock States Preparation 


Optical trap 



Figure 9 (a) A nanodiamond with a NV center is optically trapped in 

vacuum with spin-mechanical coupling enabled through a nearby magnetic 
tip and opto-mechanical coupling through a cavity around, (b) The atomic 
structure (left) and the level diagram (right) in the ground state manifold for a 
NV center in the nanodiamond. Figure adopted from Copyright (2013) 
by The American Physical Society. 

As shown in Fig. [9^, we consider a nano-diamond trapped by 
an optical tweezers in a high-Q cavity in vacuum. Therefore, 
the mechanical motion of the nanodiamond couples with the 
cavity mode. The trap is located in a place that the cavity 
mode has the maximum gradient. Near the NV center, there 
is a magnet tip, which induces a strong magnetic field gra¬ 
dient. The magnetic gradient couples the mechanical motion 
and the electron spin. There is also a microwave source to 
control the spin of the NV center inside the nanodiamond. As 
the nanodiamond is optically trapped in the vacuum, the co¬ 
herence time of the mechanical mode for its center-of-mass 
motion is long tSOll . The frequency of the optical trap can 
be quickly tuned through control of the laser intensity. This 
feature is important as we can cool the mechanical mode to 
the ground state in a strong trap and prepare large quantum 
superposition states of the nanodiamond in a weak trap by a 
quench of the trap frequency. 

Once the resonator cooled down, we can generate many 
non-classical states, such as Fock states. The NV spin is ini¬ 
tially set to the state |0), which is decoupled from the me¬ 
chanical mode during the cooling. Initialization and single 
shot detection of the NV spin have been well accomplished 
experimentally iMl. We assume that the NV center is at 
a position with zero magnetic field and a large field gra¬ 
dient. We apply a microwave drive with the Hamiltonian 
Hdrive = + Qnv,-i^^‘""10)<-11 + h.c.)/2 

and set the Rabi frequency Onv,+i = f^NV and the detuning 
A+ = COI+-CO +1 = A. With A » |OnvL we adiabatically elim¬ 
inate the level |0) and get the following effective Hamiltonian 

He = flCOmCl^nClm + + flX(cr+ + (T-){ani + ( 14 ) 

where Q = |Qnv|V4A, ct, = |+)(+| - |-){-|, = |+){-|, 

(T_ = |-)(+|, and we have defined the new basis states 
1+) = (I + 1) + I - 1))/V2, I-) = (I + 1) - I - 1»/^. In 
the limit A we set Q = and use the rotating 


wave approximation to get an effective interaction Hamilto¬ 
nian between the mechanical mode and the NV center spin, 
with the form 

Hjc = hAcr+am + h.c.. 

This represents the standard Jaynes-Cummings(J-C) cou¬ 
pling Hamiltonian. Similarly, if we set Q = the anti 

J-C Hamiltonian can be realized with 

Hajc = fiAcr+a]^ -f h.c.. 

Arbitrary Fock states and their superpositions can be pre¬ 
pared with a combination of J-C and anti J-C coupling Hamil¬ 
tonians. For example, to generate the Fock state |2)^, we 
initialize the state to |+)|0)^, turn on the J-C coupling for a 
duration t\ = 7il(2A) to get |-)|1)^, and then turn on the anti 
J-C coupling for a duration t 2 = fi/ V2 to get |-f)|2)^. The 
Fock state with arbitrary phonon number can be gener¬ 
ated by repeating the above two basic steps, and the inter¬ 
action time is ti = til for the ith step. Superpositions of 
different Fock states can also be generated. For instance, if 
we initialize the state to (co|+) + ci|-)) 0 \0)m/ V2 through a 
microwave with arbitrary coefficients cq, ci, and turn on the 
J-C coupling for a duration ti, we get the superposition state 
|-)(8)(ci|0)^-r/co|l)m)/ V2. Using the optical cavity, the Fock 
state \nm)m of mechanical mode can also be mapped to the 
corresponding Fock state of the output light field 1431 . 

The effective QND Hamiltonian for the spin-phonon cou¬ 
pling takes the form 

Hqnd = hxo-zol^am, 

with;^ = 4^A^/(4Q? - oj^) when the detuning ||Q| -6 g)^/ 2| » 
A. The Hamiltonian can be used for a quantum non¬ 

demolition measurement(QND) measurement of the phonon 
number: we prepare the NV center spin in a superposition 
state |-r) -r ^^2, and the phase 0 evolves by 0(0 = 

00 + Axrijnt, where = aham denotes the phonon number. 
Through a measurement of the phase change, one can detect 
the phonon number. 

Let us estimate the typical parameters. A large magnetic 
field gradient can be generated by moving the nanodiamond 
close to a magnetic tip. Though magnetic gradient up to 
10^ T/m has been realized ||65l[66l, here we take the gradi¬ 
ent G = 10^ T/m, which is much less than the one used in 
Sec. 12.1.21 ||29l. We get the coupling A ^ An x 52 kHz for 
a nanodiamond with the diameter d = 30 nm in an optical 
trap with a trapping frequency ojm = 2n x 0.5 MHz. The 
QND detection rate 2\x\ ~ 2n x 25 kHz with the detuning 
||Q| - 0Jm/2\ ~ 5A. The NV electron spin dephasing time over 
1.8 ms has been observed at room temperature ll34l , which is 
long compared with the Fock state preparation time 1/A and 
the detection time 1/ (2\x\). 

3.2 Macrocopic quantum superpositions 

To prepare spatial quantum superposition state, we need to 
generate quantum superposition of the nanodiamond at dis- 
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tinct locations. To detect the superposition state, we need to 
do interference experiment, either with matter wave ll50l or 
with the spin ll52ll . 


3.2.1 Matter wave interference 


Without the microwave driving, the spin-mechanical cou¬ 
pling Hamiltonian takes the form 


H = fi(jJma]^am + JiASzicim + ajj. (15) 


The mechanical mode is initialized to the vacuum state |0)^ 
(or a Fock state in a strong trap with the trapping fre¬ 

quency ojmo and the NV center spin is prepared in the state |0). 
Although the ground state cooling is most effective in a strong 
trap, to generate large spatial separation of the wave packets 
it is better to first lower the trap frequency by tuning the laser 
intensity for the optical trap. While it is possible to lower 
the trap frequency through an adiabatic sweep to keep the 
phonon state unchanged, a more effective way is to use a non- 
adiabatic state-preserving sweep ll67]| . which allows arbitrar¬ 
ily short sweeping time. We denote as the mechanical 

state in the lower frequency ojmi . We then apply an impulsive 
microwave pulse to suddenly change the NV spin to the state 
(|-rl)-r|-l))/V2 and simultaneously decrease the trap fre¬ 
quency to ojmi ^ ojmi . The evolution of the system state under 
the Hamiltonian (fTSl) then automatically split the wave packet 
for the center-of-mass motion of the nanodiamond. The split¬ 
ting attains the maximum at time T 2 I 2 = where the 

maximum distance of the two wave packets in the superpo¬ 
sition state is Dm = ^Aa 2 loJm 2 = AgsiaBGI(m(o^ 2 )^ where 
a 2 = AlfillmcDmi^ At this moment, the system state is 


I'Ps) = (I + l>|£>m/2)„„ + I - 1)1 - DJ2),JI ^l2, (16) 


where | ± D„l2}n^ = \n„}i is the dis- 

placed Fock state (or coherent states when nm = 0). This is 
just the entangled spatial superposition state. 

The maximum distance Dm is plotted in Fig. [TOh versus 
trap frequency, magnetic field gradient, and diameter d of the 
nanodiamond, and superposition states with separation Dm 
comparable to or larger than the diameter d is achievable un¬ 
der realistic experimental conditions. 



Position(nm) G(T/m) 


Figure 10 (a) Spatial interference patterns for a 30 nm nano-diamond af¬ 

ter 10 ms of free expansion. The nano-diamond is initially prepared in the 
vacuum state 10)^ or the 1-phonon state |l>m of a 20 kHz trap. The magnetic 
gradient is 3 x lO^T/m. Before the trap is turned off, the center of mass of the 
nano-diamond is prepared in |j^+>o. (b) Maximum spatial separation Dm as 
a function of the magnetic gradient G when the trapping frequency is 1 kHz. 
Macroscopic superposition states with separation larger than the size of the 
particle can be achieved with a moderate magnetic gradient. Figure adopted 
from Eo). Copyright (2013) by The American Physical Society. 

To transform the entangled cat state I'Ty) to the standard 
cat state |(A±>„„, = (|£)„/2)„,, ± | - £)„/2)„„)/ yl2, we need 
to apply a disentangling operation to conditionally flip the 
NV spin using displacement of the diamond as the control 
qubit. This can be achieved as different displacements of the 
wavepacket induce relative energy shifts of the spin levels 
due to the applied magnetic field gradient ll50ll . To detect spa¬ 
tial superposition state, we can turn off the optical trap and 
let the spatial wave function freely evolve for some time t. 
The split wave packets will interference just like the Young’s 
double slit experiment. The period of the interference pat¬ 
tern is Az = Anfit/irnDm). As an estimation of typical pa¬ 
rameters, we take ojmi = comi = An xlQ kHz, J = 30 nm, 
and magnetic field gradient 3 x 10^ T/m. The spin-phonon 
coupling rate A ^ In x 11 kHz and the maximum distance 
Dm - 31(22. The preparing time of superposition state is 
about 25 yus, which is much less than the coherence time of 
the NV spin. For the time of flight measurement after turn¬ 
off of the trap, we see the interference pattern with a period 
of 47 nm after t = 10 ms, as shown in Fig. [TOl which is large 
enough to be spatially resolved unma. 

3.2.2 Ramsey interference 

The challenging of matter wave interference proposal is that 
the mass and velocity various for an ensemble of nano¬ 
diamond II 52 I . In order to overcome this challenging. Ref. 
II 52 I proposed to used Ramsay interferometry to observe the 
interference. With Ramsay interence, we can verify the mat¬ 
ter wave superpositions, as well as the quantum contextuality 
for macroscopic quantum systems ll5^ . The basic setup of 
their poposal is the same as Fig. |9l Here we assume that an 
angle 6 between the vertical and the z axis of the system. The 
Hamiltonian of the system is 

H = fiDS^ + hojral^am - 2fi(ASz - AA)(am + < 2 ^), (17) 
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where D = 2.87 GHz is the zero field splitting, and Ad = 
^mg cos 6 . Let’s denote ii = 2(5^^ - lAX)lfi(jJr- The 

initial state of resonator is a coherent state |yS). The spin 
dependent time evolution of the state is \l3{t,Sz))\Sz), where 
Sz = 0, +1. We get that 

Pit, SP) = ^ 


A striking feature of this evolution is that the oscillator re¬ 
turns to its original coherent state p for any p and at time 
t{) = InjcOr^ 

As an example, we consider a initial separable state |yS)|5'^). 
By using a strong microwave pulse H^w = + 1)(0| + I “ 

1)(0| -r h.c., we can change the state to |'L(0)) = P)(\ + 1) + 
I - 1))/ V2. The state after an oscillation period to = InjcOr 
becomes 

J + 1) + _ 1). 

m/o)> = —A__J—(18) 

where we have Acpg = \6AAAtol(fi^(ji>r)- We apply microwave 
Hamiltonian H^w again to transfer the spin population to |0). 
After time tp = 7r/(2 V2Q), we get 


PO = COS2(^), 


which can be used to measure the relative phase Acpg. The 
detection of A(pg is an evidence for the matter wave superpo¬ 
sition state. Besides, as the phase A(pg is not dependent to the 
initial state p, it is possible to be observed when initial state 
is thermal. 


3.3 Applications as mass spectrometer 

As mentioned above, one of the most attractive features of the 
system of optically levitated particles is its high mechanical 
quality factor. The trapped particle can coherently oscillator 
many times (e.g., 10^^ times) before its mechanical energy 
is damping out. When this kind of almost-perfect mechani¬ 
cal oscillator couples to a long-live quantum object, i.e. NV 
center electron spin, fantastic quantum phenomena and novel 
application can be expected. 

In Ref. 1541, we studied the quantum dynamics of a cou¬ 
pled system of a quantum oscillator and a single spin, where 
the spin is under dynamical decoupling control. Dynamical 
decoupling control, as an extended version of spin echo tech¬ 
nique, has been proved to be an efficient method to prolong 
the spin coherence time. In Ref. we predicted that, under 
dynamical decoupling control and coupled to a high-quality 
mechanical oscillator, the spin coherence will exhibit a se¬ 
ries of periodic peaks, forming a comb structure (called time- 
comb in Ref. US). The peak width is calculated as 


A 


q* - 


2 V 2 

r?' 


To 

NA V2«th + 1 


(19) 




time (irs) time (irs) 


Figure 11 (a) The time-comb structure of qubit coherence under 100-pulse 
CPMG control. For ojrt » 1, the comb period is synchronized with the os¬ 
cillator period Tq = Injajr- (b) Close-up of the coherence peaks. A missing 
peak 2ii (jjyt = Nn is indicated by the blue dashed line. The peak width is 
decreasing when getting close to the missing one. (c) Close-up of the nar¬ 
rowest coherence peak, which is centered at tq* with width Aq* (see text). 
The parameters used in this figure are oscillator frequency 100 kHz, cou¬ 
pling strengt A = O.OOOluj^, temperature T = 10 K and 100-pulse CPMG 
control. 


Two features are notable in the above equation. First, the 
peak width is inversely proportional to the square root of the 
thermal occupation number ^th of the mechanical oscillator, 
which implies narrower peaks in higher temperature. Second, 
the peak width is inversely proportional to the control pulse 
number N. Since the control pulse number N is related to the 
total coherent evolution time, and the later can be regarded 
to a kind of resource in high-precision sensing, we expected 
that improving the dynamical decoupling control pulse num¬ 
ber can improve the measurement sensitivity. 

Based on the time-comb structure, we propose an 
ultra-sensitive mass spectrometer using the coupled qubit- 
oscillator system. Since the peak position in the time axis is 
synchronized with oscillator period, we propose to infer the 
oscillator period (and its mass assuming a constant spring co¬ 
efficient) by measuring the peak position. Obviously accord¬ 
ing to the error propagation formula, the narrower the peak 
is, the better sensitivity we will have. The two features men¬ 
tioned above make the proposed mass spectrometer different 
to those based on traditional measurement principles. As the 
peak will be narrower in higher temperature, we propose that 
the mass spectrometer will have a counterintuitive temper¬ 
ature dependence of the sensitivity (i.e., higher temperature 
helps improve the sensitivity). This result makes the room- 
temperature application of ultra-sensitive mass spectrometer 
possible. The second feature, as we showed in Ref. ll54l . im¬ 
plies an improved scaling of the measurement sensitivity to 
the measurement resource. The scaling relation goes beyond 
the shot-noise scaling and Heisenberg scaling in the quan¬ 
tum metrology theory. The underlying physics which brings 
about this improved scaling is under further studying. 
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3.4 Experimental progresses 





Figure 12 A photograph of a ~ 100 diameter nanodiamond which is lev¬ 
itated by a optical tweezer (1064 nm). The green color is from an optical 
pump 532 nm laser. The trapping lens is on the left, and the lens on the 
right is used to collimate the light exiting the trap. Figure adopted from ES). 
Image credit: J. Adam Fenster, University of Rochester. 

Dielectric nanoparticles have been trapped and cooled opti¬ 
cally, as shown in Ref. 114111421 l44ll . Following the similar 
method, trapping nanodiamond in atmosphere pressure has 
been realized experimentally in 2013 ED. The experimental 
setup is shown in Fig. [121 They found the evidence of NV 
photoluminescence from a nanodiamond in a free-space op¬ 
tical dipole trap. The photoluminescence rates are decreased 
with increasing trap laser power. The neutral charge state 
(NV^) can be suppressed by continuous-wave trap. 

The charged nanodiamond was also levitated in a linear 
quadrupole ion trap in air under ambient conditions |[68l . An 
electric trap has no trapping laser that quenches fluorescence 
emission. The surface charge has effects on the emission rates 
of fluorescence spectrum of NV center. However, there is no 
difference in the spectral properties was found in the experi¬ 
ment. Up to now, there is no report on trapping nanodiamond 
in vacuum. This is the first problem that should be overcome 
in future experiment. 

4 Conclusion and Outlook 

In this paper we gave a brief review on the held of 
(opto)nano-mechanical resonator coupling with NV centers 
electrons spins. We discussed magnetic and strain induced 
phonon-spin coupling mechanisms. We discussed how to 
cool the resonator to the ground state and how to manipu¬ 
late the phonon non-classical states. We also discuss how to 
realize quantum information processing in this system. We 
reviewed the progresses in optically trapped nano-diamond 
with building in NV centers. We discussed how to observe 
macroscopic quantum interference in the system. We also 
discussed how to used the system as an ultra-sensitive mass 


spectrometer in room temperature. Finally, we briefly sum¬ 
marized the experimental progresses. 

Here we discuss some new directions and ideas in the held. 
Up to now, we only consider the first order magnetic gradient 
induced coupling. It would be interest to study the second or¬ 
der gradient induced spin-resonator coupling. This coupling 
mechanism may have some new applications and features. 
For example, it may be used to entangle NV centers with dif¬ 
ferent transition frequencies. It may also be used for cooling 
mechanical resonator which is out of Lamb-Dicke regime. 

Another new direction is to study the coupling between 
rotation degree of freedom of trapped nano-diamond to the 
internal NV centers electron spins. Geometric phase, even 
non-Abelian geometric phase may be induced by mechanical 
rotation (6373. The rotation may also be cooled and manip¬ 
ulated by NV centers. 

It would be interesting to investigate how to cool down the 
internal temperature of the diamond which hosts the NV cen¬ 
ters. As we know, the decay rate of the NV centers electron 
spins is in verse proportional to the temperature. Only under 
the low temperature T < 10 K, the decay rate can be de¬ 
creased to around Hz. By considering the coupling between 
NV center electron spins and the phonons in the diamond, 
it would be actively cooled down the temperature of a nan¬ 
odiamond trapped with an optical tweezers ||7Tl|72l . In this 
way, we may extend the lifetime of NV centers electron spins 
around second even at room temperature environment. 

Before conclusion, we note that the proposals and schemes 
discussed in this review should also be applied in other simi¬ 
lar systems, such as hybrid mechanical resonator with quan¬ 
tum dots |[73ll . ^8Si nanoparticles with donor spins |[74ll . or 
nanocrystals doped with rare-earth ions liTSl . 
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